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ABSTRACT 

The Naval Post Graduate School has obtained software which can be used for the 
thermal analysis of an electronic component. This software includes two versions, one for 
steady state analysis and the other for transient analysis. Each version consists of two 
programs, a model builder and a thermal analyzer. This thesis describes a user friendly, 
menu driven specific model builder which may be used to rapidly generate a thermal model, 
containing up to 750 nodes, for a microcircuit die. This model builder is an improvement 
over the existing model builder in that the only pertinent input will be the physical 
dimensions and heat transfer data. The main feature of this model builder is the 


accomodation of the marked variation of silicon thermal conductivity with temperature. 
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l. INTRODUCTION 

Microcircuit chips have become commonplace items in all types of electronic 
equipment, both military and commercial. The capabilities of these chips continue to 
increase as improvements in design and manufacturing are incorporated. However, the 
design process becomes more difficult as the ability to place more and more circuits on 
these chips is realized. One area of immediate concern is in thermal control. Chips 
containing more than one million components in a 60 mm? area are known to exist [Ref. 
1). Heat dissipation from components contained in the chip as well as from any outside 
heat sources needs to be thoroughly analyzed. The reasons for conducting a thermal 
analysis on a microcircuit structure concern the impact of high temperature on reliability, 
bias stabilization, manufacturing costs, and thermal runaway. 

The reliability of a microcircuit is directly related to its junction temperature. It is now 
common knowledge that failure rate is proportional to junction temperature and designers 
plan for the microcircuit to operate within a specified temperature range. If this range is 
exceeded, the operating capability of the device may be severely reduced or it may fail 
completely. Because the failure rate increases exponentially as the temperature increases, 
the effects of high junction temperature can be drastic. 

It is necessary to bias an electronic circuit in order to obtain a stable and predictable 
operating point. As the junction temperature changes, the impact on this operating point 


will be directly impacted. Consider Figure 1 and observe that by Kirchhoff’s voltage law, 


-Vec + Re Ie + Vee + Re I, = 0 (1) 
so that the collector to emitter voltage will be 
Vee = Vee - Ic (Re + Re) (2) 
1 














Voc = power supply voltage to collector 
B = base 

C =collector 

E =emitter 

Voce = voltage across collector and emitter 
R, = emitter resistance 

R, = collector resistance 

I, = emitter current 

I , = collector current 


Figure 1 

















It is a fact that, the collector current, I, is directly proportional to the junction temperature. 
As the junction temperature increases, the voltage drop across R, and R, increases and as a 
result, Ve is quite likely to decrease to an intolerably low level. 

During manufacturing, the chip, at some point, must undergo a soldering process in 
which connections are made from chip to chip or from one component to another. A 
particular soldering compound must be selected which will accomodate the expected 
operating temperature range without failure. The soldering compound cost is based on it’s 
composition and, as general rule of thumb, the higher the melting point, the higher the cost. 
An increased cost incurred by a manufacturer for a high temperature solder may be very 
minimal on a per unit basis. But any manufacturer dealing with microcircuit chips may be 
handling millions of them and these costs could conceivably reach the million dollar mark 
which is probably not recoverable. 

Hot spots anywhere in the microcircuit chip could result in thermal runaway. Because 
the temperature increases as the collector current increases, a hot spot near the transistor 
will cause an increase in the collector current. This then causes an increase in power 
dissipation (P = I, Ve) which in turn, causes an increase in temperature. This increase in 
temperature causes an additional increment of collector current to flow and leads to yet 
another increase in power dissipation. This regenerative process is known as thermal 
runaway and as the cycle continues it could lead to drastic and complete device destruction. 
(Ref. 1] 

The Naval Post Graduate School has obtained software which can be used for the 
thermal analysis of an electronic component. This software includes two versions, one for 
steady state analysis and the other for transient analysis. Each version consists of two 
programs, a model builder and a thermal analyzer. The model builder is designed to 


provide the user with the means to model a physical configuration in order to provide the 








thermal analyzer with an input file. This file is taken by the thermal analyzer as its input 


and an output file containing a summary of the temperatures within the configuration will 
be generated. 

The generation of the thermal model is the biggest hurdle in the actual analysis and the 
analyst is powerless in the absence of a model builder program. However, the model 
building process begins with a drawing of the configuration and a subdivision of it into 
small but finite subvolumes. Each subvolume is presumed isothermal and the centers of 
each subvolume are then representative of the entire subvolume. These centers (centroids) 
are called node points or merely nodes and thermal analysis involves the writing of n node 
equations in nm unknown temperatures where the nodes are connected by thermal 
conductances. 

The very laborious and time consuming process of writing the node equations can be 
avoided through the use of a model builder which will write them automatically. Yet, there 
are drawbacks in the use of a general model builder in that, while it can be used for a 
specific geometry and set of material conditions, many hours can be spent on input in front 
of a computer terminal and, before using any model builder, careful consideration must be 
given to whether or not the configuration under analysis will be sufficiently well 
represented by the model generated. 

This thesis describes a specific model builder which may be used to rapidly generate a 
thermal model containing up to 750 nodes for a microcircuit die. The main feature of this 
model builder is the accomodation of the marked variation of silicon thermal conductivity 
with temperature. 

The specific model builder has additional features. These will be described in this 


thesis and are as follows: 


1) The capability of working in SI or English units. 

2) The user has five aspect ratios from which to choose. 

3) The provision of up to 750 nodes for the die model. 

4) The ability to input heat dissipation using several methods. 

5) The automatic calculation of conductance values based on user input 
dimensions. 

6) The automatic correction for variations in the thermal conductivity values of 
silicon with temperature and location. 


7) The provision of automatic linearization of radiation. 














il. A PROBLEM IN THERMAL ANALYSIS 

The provision of a solution to a simple problem in heat transfer analysis will show that 
a thermal analyzer is required for serious, large scale, thermal analysis. 

Conduction is the heat flow mechanism whereby heat is transfered from one part of a 
homogeneous medium under the influence of a temperature gradient without appreciable 
displacement of the particles that compose the medium. It was Fourier who proposed that 
the heat flow is directly proportional to the area of the heat flow path and the temperature 
gradient along the path, 

qe ag B) 
Insertion of a proportionality constant yields: 

=-kA a (4) 

where 

A = area of the heat flow path 

k = thermal conductivity of the material 

a = change in temperature per unit path length (temperature gradient) 
where the minus sign assures a positive heat flow in the presence of the required negative 
temperature gradient. [Ref. 2] 

Equation (4) serves to define the thermal conductivity of the material 
k=- “at W/m-°K (5) 
dx 

and it must be noted that the simplicity of this development quickly changes when k is a 
function of temperature, as in the case of silicon. 

Equation (4) is a simple differential equation which can be solved using the 


boundary conditions in Figure 2. Separation of the variables gives: 






































- pk dx =dT 

and the i ecg perwern the limits: 
“h Jax = = jar 

gives: 

(x2 - x1) =T2-T; 
With L = xg - x), it is seen that: 

q=“Aq,-T» (6) 
which is called Fourier’s law. 


Consider, Ohm’s law which governs the flow of current in an electrical resistor: 


AV 
= ) 


and note the similiarity between the forms of equations (6) and (7). An analogy may be 


proposed where the analogous quantities are: 


qe! (8a) 
AT > AV (8b) 
& eR (8c) 


In particular, equation (8c), defines the thermal resistance for the conduction process and 
one observes that, for the slab in Figure 2, if T, and T2 are given, then q can be determined 
if the thermal resistance of the slab is known. On the other hand, if T,; and q are known, 
T2 can be calculated using R = rx 
T2=1T,-Rq 

Convection is a fluid flow process and, if it applies in heat transfer, it is the transfer of 
heat from or to a confining surface by a flowing fluid. The fluid flow may be induced by 
buoyancy or density gradients within the fluid and this phenomenon is called natural or free 
convection. If the fluid flow is induced by a fan, a pump, or a blower, it is termed forced 


convection. [Ref. 2] 








Heat flow by convection is directly proportional to the temperature difference between 


the confining surface and fluid and the surface area over which the process takes place: 

qx AAT (9) 
and a proportionality factor may be inserted: 

q=hAAT (10) 
Equation (10) is the Newton law of cooling and h is called the coefficient of heat transfer 


defined by: 


h =—4— W/mK (11) 
A AT 
Notice that the concept of thermal resistance can be extended to the convective case 
where: 
~ =R (12) 


If convection exists on both faces of the slab in Figure 2, the configuration of Figure 3 


exists and the electric analog is shown in Figure 4.Here the total resistance is: 


AT rT, -T 
ies geen) 1 Se 
R (13) 


Thus far, no difficulty has been encountered in a calculation procedure. The picture 
changes dramatically when radiation is considered. 
Radiation is electromagnetic in nature and is the means for heat communication in a 


vacuum. It is beyond the scope of this thesis to derive the general radiation equation: 














Figure 4 








q=0F, F, A (1,4 - T,4) (14) 
where 
6 is the Stefan-Boltzmann constant = 5.67(10-8) W/m2-K4 
Fg is a shape factor that accounts for the arrangement of the radiating source and 
the absorbing receiver 
F, is the emissivity factor that accounts for the ability of the surfaces to absorb or 
emit radiation 
T; is the temperature of the source, Kelvin 
T, is the temperature of the receiver, Kelvin 
A is the surface area, m2 
In order to put the radiation on a q = K AT basis where K is the thermal conductance, 
consider that: 
T,4 - T,4 = (T,? + T,) (T,? - T;*) 
and: 
T,4 - T,4 = (T,? + T,?) (T; + T, (Ts - Ty) 
Equation (14) may then be written as: 
q=6 Fag Fe A (1,2 + T,*) (T; + Ty (T; - Ty) (15) 
and a radiative heat transfer coefficient may be defined: 
q =h, A (T; - T;) (16) 
where 
h, = 0 Fag Fg (T,2+T,*) (T;+Ty) (17) 
A thermal resistance in radiation can then be proposed as: 
R=iK (18) 


and, as shown in Figure 5, this resistance is in parallel with the convective resistance. 
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q kA 
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ee ae Ree 
h;, A h A 
T, 


2 2 

h = OF, F. (T, +T, CT, +1 ) 
2.2 

h, = OFA F. (T, +T, ) (T, + T, ) 


Figure 5 
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A solution procedure for the thermal network in Figure 5 must be iterative because of 


the form of h, given by equation (17). This is obviously no longer a simple problem and, 


at this point, one must turn to the computer for help with more complicated configurations. 
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itl. FINITE DIFFERENCES AND NODAL ANALYSIS 
It is the purpose of this section to provide the assurance that the method of modelling 
will yield accurate results and to show how the n node equations in the n unknown 
temperatures are formulated. 
Equation (4) is but a one dimensional simplification of the more general equation of 


heat conduction [Ref. 3]: 


d 7. aT, . Of. ary 28 yar aT 
x (k5x) +3, (k3y) ++ (kS) + qi = pc 5 (19) 


where 

T = temperature, °F or °C 

t = time, hrs. 

qj = rate of internal heat generation, Btu/hr or Watts 

k = thermal conductivity, Btu/hr-ft-°F or Watts/m-°C 

P = specific weight, lbm/ft3 or kg/m3 

x, y, Z = coordinate system distances, ft or m 

c = specific heat, Btu/lbm-°F or Watt-hr/kg-°C 
We note that if we wish to formulate a problem that deals in discrete Ax, Ay, Az, and 
AT’s, we must be able to accurately approximate the first and second derivatives. 

In general we must be able to determine the slope of a line tangent to the point of 

interest on the plot of the temperature as a function of x, y, or z. (Figure 6). Using x as the 
spacial coordinate, we observe that the slope of the line (m), if it exists, will be: 


we Ake Sot Ax) - F(Xo) (20) 


Ax 
The term, [f(x,+Ax) - f(x,)]/Ax, is called the difference quotient and is the ratio of the 
change in the value of the function at x, and x,+Ax to the change in x . The limit of the 


difference quotient is called the derivative of the function at xq: 
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T=F(x) %+Ax 


Figure 6 
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lim ot Ax) - f(Xo 
F) = Ae F(Kot ~ F(Xo) (21) 


In examining Figure 6 we can derive the derivative for this particular curve at the point x 
by letting f(x) = T(x) and noting that: 


f(Xo+ Ax) = T(Xp+Ax) 


f(Xo) = T(Xo) 
i ° - T(Xo 
f’ (Xo) a one _T(Xo+Ax) - T(Xo) (22) 


Ax 
and when Ax — 0 we observe that the expression for the first derivative becomes: 


iy) (23) 
Ax 
and we note that: 
AT 
4 =—— 24 
f'(%o) ie (24) 


as long as Ax is small as required by the limiting process. [Ref. 4] 
The second derivative, using the defintion of a derivative, will be the ratio of the 
change in the value of the first derivative to the change in Ax. Thus, the second derivative 


in Figure 6 becomes: 
f'(X) = T(X9+Ax) - T(Xo) 
Ax 


T(X9+2Ax) - T(X9+Ax) 


Y ot = 
Ff’ (Xt Ax) a 


im _f’(Xot Ax) - f'(Ko 
fx) = gin, Feats) FH) 


Ax 
F°(xo) = Ttot2An) - TexotAx)  T(xotAx) - T(%e) 
Ax Ax 
Ax 
#°Ciq) = Tot ZAW) = Tot AX) + To) 2 


Ax2 
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As in the case of the first derivative, the second derivative given by equation (25) is a 


valid approximation as long as Ax is small. The use of a Taylor series expansion will 
confirm the validity of equation (25). Taylor’s theorem states that a function can be 


approximated by a polynomial of the form: 
(n) 
f(x) = f(a) + f’(a)(x-a) + LO. aj? +...+ ayn (26) 


where the polynomial is for the function expanded about x = a [Ref. 4]. Using our 


example the Taylor series for the function T(x) at T(x9+2Ax) is: 
2 
STxotAn) Aug 1 d2T(x+Ax) 


T(x,+2Ax) = T(xg+Ax) + 5 dx? 


Ax2 + 


1 d3T(x,.+Ax) 4 La dT (Xo+Ax) 
6 dx3 AX? +... + 7 dx" (27) 
and for the function T(x) at T(x,): 
7 AT(Xo+AX) , 1 PT (Xo+Ax) 4 
T(Xo) = T(Xo+Ax) - Gy — Ax + a2 ~~ 7) 
1 d3T(x9+Ax) 3 mes dT (XotAx) 
6 a Ax3 +. = al axa 4* (28) 


If we keep Ax very small then the terms smaller than Ax? become negligble. If equations 


(27) and (28) are added we can obtain an expression for the second derivative: 


2! 
T(Xo+2Ax) + T(Xo) = 2T(Xo+Ax) + oe Ax? (29) 


d2T(xo+Ax) _ P(Xo+2Ax) - 2T(X9+ Ax) + T(Xo) 
dx? Ax2 
which confirms our earlier result. 


(30) 


Thus, as long as Ax is kept very small, the following are valid representations and can 


be used in the analysis process: 
dT _ T(X+Ax) - T(Xo) _ AT 


= = (31) 
dx ~ Ax Ax 

d2T T(x 9+2Ax) - 2T(x,+Ax) + T(x) 32 
qn ge ry 
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As stated in section I, the configuration to be analyzed can be divided into small but 
finite volumes which are considered to be isothermal. The geometric center of each 
subvolume is designated as a node and the node can be representative of the entire 
subvolume. These nodes may be connected to each other through thermal resistances and 
the nodal analysis for the node temperatures lends itself to solution on the computer. 

The first law of thermodynamics (the law of the conservation of energy) states that 
energy cannot be created nor destroyed, but can be transformed from one form to another. 
One may therefore form an energy balance on the typical node shown in Figure 7. For ease 
of explanation, the environment will be considered to be a single node, node 100 (Figure 
8). Examination of node 5, shows that it is connected to nodes 2, 4, 6, 8, 14, and 100. If 
any energy is directly applied to or removed from this node, we would also include it in the 
equation. The node equation for node 5, with rate of heat input (qj), becomes: 

q2 + G4 + Go + G8 + G14 + Gioo - Gi = 0 (33) 
where each one of the q’s with a numerical subscript represents the rate of heat flow from 
node 5 to the node indicated by the subscript. 


Equation (33) is a continuity relationship and it does not show the items of interest, the 


node temperatures. But q = at and the terms in equation (33) can be given by: 
kA k Ax A 
q2 = AT =———— (Ts- T2) (34a) 
Ay 
kA k Ay Az 
qa= 9 AT =—— = (T5- Ta) (34b) 
Ax 
kA k Ay Az 
qo=7 AT =U (T5- Te) (34c) 
Ax 
kA k Ax Az 
qg=7 AT aie one (Ts - Tg) (34d) 


19 














T1009 is on the top face of node 5 and, 


Ax, Ay, Az = the distance between nodes in the x, y, and z directions respectively 
If we assume that Ax = Ay = Az = 1 the node equation becomes: 
K(T5-T2)+K(Ts-T4)+K(T5-T6)+K(T5-Tg)+K(T5-T14)+2K(T5-T100) = qi (35) 
or 
-Tz - Tq + 7Ts - Tg - Te - T14 =f + 2T 300 (36) 
where the terms on the right hand side of equations (35) and (36) should be known values. 
A node equation can be written for each node in the configuration and in Figure 8 
which contains 18 nodes, 18 node equations would be required to obtain the 18 unknown 
node temperatures. To make maximum use of the thermal analysis software, the generation 
of 750 nodes with 750 node equations similar to the one just derived would be required. 
There are a number of methods that can be utilized to solve these equations. It is not our 
purpose to solve these equations at present but to input them into a computer program 


which will solve them for us. 
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IV. MICCKT AND THE THERMAL ANALYZER 
The model builder program which will become part of the thermal analyzer software 
package is MICCKT. MICCKT is for the specific microcircuit configuration shown in 
Figure 9, specifically, the microcircuit is covered by a lid and supported by a base. The 
gap between the chip and the lid is assumed to be air filled. Once the physical properties of 
the microcircuit are entered MICCKT will generate an output file in the correct format for 


use in the thermal analysis software. 


A. MICROCIRCUIT MODELLING PROCEDURE 

Recalling the idea that Ax, Ay, and Az must be kept small to enable the use of the finite 
difference approach, we divide the microcircuit into subvolumes. The model builder will 
divide the chip into four layers with the top and second layers containing four times as 
many subvolumes as the bottom layers. (Figure 9) The top layer, containing the circuitry, 
is preset at 0.002 inches (0.00508 cm) and the thicknesses of each of the next three layers 
are one third the overall thickness less 0.002 inches. The reason for this is that heat 
injection will occur in the top layer and accuracy is assisted when the spatial increment is 
smaller. In accordance with the procedures already discussed, the node will be placed in 
the geometric center of the subvolume and will be representative of the physical properties 
of that subvolume. 

The microcircuit structure is made of pure silicon. This presents a significant 
challenge as the thermal conductivity of silicon varies widely with temperature. (Figure 10) 
In order to determine the internodal conductances it is necessary to modify this curve for 
use in the thermal analyzer. The thermal conductance is function of the area of the path of 


transfer, the path length, and the thermal conductivity. Because the area and path length are 
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Thermal Conductivity, Ke, 


Wi(em —°C) 


IMPURITY 


no added 
impurity 
As 
Au 


unknown 
B 
8 
unknown 


p-Type 





Temperature (K) 


(MPURITY 
CONTENT 
(em~3) 


CROSS SECTION 
(mm X mm) 


2.15 x 1.98 


unknown 
2.36 x 2.21 
1.75 x 1.5 
3x3 
3x3 
unknown 
unknown 





Figure 10 [From Ref. 5] 
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established by the internodal distances, the thermal conductivity at a given temperature is 
simply multiplied by a scalar value (A/L). It is the selection of the thermal conductivity for 
the desired temperature that poses a significant problem. 

One of the features of the thermal analyzer is its ability to generate a limited facsimile 
of a temperature dependent curve based on user inputs. This temperature dependent curve 
may be specified using eight straight line segments between nine data points. Each data 
point represents a temperature and the corresponding conductance for the physical and 
thermal characteristics of the chip. Because the physical dimensions of each subvolume 
may change from location to location, several curves may be required to completely analyze 
the configuration. 

It was determined that eight curves representing the conduction mode of transfer 
within the structure were required for the model. The top layer required two curves, 
designated curves 1101 and 1102, to represent the north-south and east-west node 
connections and a third curve, 1103, to connect the nodes in layer one and layer two. The 
second, third, and fourth layers use the same curves, 1104 and 1105, for the north-south 
and east-west connections. The Ax and Ay distances in the bottom two layers are twice 
that in the second layer but because either Ax or Ay will represent the path length, the 
curves are identical. Curve 1106 is used to connect layer two to layer three and curve 1107 
to connect layer three to layer four. Curve 1108 is used to connect the fourth layer nodes to 
the heat sink. 

The generation of each curve requires the establishment of the thermal conductivity at 
specific temperatures. Since there is a nine point limitation, the temperatures were selected 


to correspond to the expected operating range of the chip. The points selected were: 
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°C Thermal conductivity W/cm-°C 


0.0 1.68 
25.0 1.49 
60.0 1.27 
85.0 1.14 

120.0 1.01 
160.0 0.88 
205.0 0.77 
255.0 0.66 
319.0 0.55 


These points were selected because they yield a piecewise continuous curve that closely 
resembles the original curve. 

The thermal interaction between the top of the chip and the lid of the container must be 
addressed. The modes of heat transfer involved are, conduction and radiation across the air 
layer. The linearization of radiation by the analyzer leads to excessive running time. Thus, 
it would be beneficial to eliminate the non linear effect of the radiation if at all possible. A 
comparision of the heat transferred by conduction to that transferred by radiation will 
provide additional information to assist in justifying an alternate approach to the inclusion 
of radiation. 


The rate of heat per unit area for conduction through the air gap is: 


Y bond = 4 (Ts - Tr) G7) 
and for radiation: 
9 bad = 6 Fa Fe (Ts? + 7,2) (Ts + Ty) (Ts - Ty) (38) 


and a comparison of these for a typical case leads to a means for the determination of the 


relative importance of each in the heat transfer process. For example, using the values: 








T; = temperature of the chip = 400 K 
T, = temperature of the lid = 300 K 
k = thermal conductivity of air at 350 K = 2.97(104) W/cm-K 
L = path length (air gap distance) = 0.254 cm 
F, = shape arrangement = 1.0 
F, = emissivity factor, for parallel plates = 
where 
€c = emissivity of the chip = 0.7 
EL = emissivity of the lid = 0.8 
the rate of heat per unit area for conduction becomes: 
9 bond = 1.17(10-) W/em? 
and for radiation: 
J haa = 5.9(10-2) Wem? 
For the selected conditions, the value for conduction transfer is approximately twice that for 
the radiation transfer. The difference between the heat transferred through the silicon and 
the heat transferred across the air layer will be on the order of several magnitudes. In this 
example, for a chip 0.05 cm thick at 400 K, the heat transfer into the second layer of the 
chip would be 109.75 W/cm2. Therefore the vast majority of heat transfer occurs 
downward into the chip. Due to it’s poor heat transfer characteristics, air essentially 
performs as an insulator. Because radiation is even less effective as a heat transfer 
mechanism the impact on the final solution is minimal. 
Although radiation effects are a small percentage of the total heat transferred, these 
effects will not be neglected in the model builder. It will be incorporated into the thermal 
conductivity of air through the use of a correction factor. The correction factor will be a 


ratio of the magnitudes of the heat per unit area for radiation to the heat per unit area of 


28 


conduction across the air layer. This correction factor can be used to derive a thermal 
conductance for heat flow across the air gap at a given temperature difference that involves 
both modes of heat transfer: 

= KA Ratio) 

Bearing in mind that the thermal transfer across the air gap by conduction and radiation 
is significantly smaller than that through silicon, and the desire to limit the running time of 
the thermal analyzer, we choose to linearize the radiative transfer mechanism in this 
manner. Conduction, with the correction factor added, will represent the entire heat 
transfer mechanism between the chip and the lid. 

Because the thermal conductivity of air is temperature dependent, an additional curve 
must be incorporated into the model builder. This curve will be designated as curve 1109 
and will be implemented in the same fashion as the silicon conductance curves. The 


correction factor will, of course, automatically be entered into this curve. 


B. FEATURES OF THE THERMAL ANALYZER 

The purpose of the model builder is to provide input data to the thermal analyzer. This 
data must be organized in a format which is consistent with the format of the thermal 
analyzer. The input file will consist of five data lines and as many as seven data sets. 
Detailed information regarding the input data file may be obtained in the user’s manual. In 
this section we describe features of specific interest as they pertain to the model builder. 
This will help the reader gain an understanding of the features of the thermal analyzer . 

Line 1 is a title line which may contain as many as 79 alphanumeric characters. This 
title will be used to identify the output files. 

Line 2 is the problem data line and will describe inputs pertaining to the problem at 
hand. Of the items in this line, the user will control the number of nodes under 


consideration and the units used. Other items identify the capabilities of the thermal 
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analyzer to handle situations involving heaters, modes of heat transfer with unique 


exponents, nodes which receive secondary heat inputs, nodes controlling fast heat, and the 
number of temperature dependent heat input curves. These entries are not applicable to this 
model builder. The number of constant temperature values and number of temperature 
dependent curves are applicable and will be set by the model builder. 

Line 3 is beyond the user’s control and is part of the thermal analysis software; all 
entries in this line are automatic. 

Line 4 is the problem capability line. This line defines the maximum values that may 
be used for each item in line 2. Of particular interest are the maximum number of nodes to 
be considered, 750, the maximum number of temperature dependent curves, 15, and the 
maximum number of constant temperatures, 50. The significance of entering the 750 node 
maximum becomes clearer when we look at the node assignments to the constant 
temperature inputs. The first constant temperature will be assigned to node 751 and 
subsequent constant temperatures will be assigned in order up to node 800. 

Line 5 items place certain controls on the computational processes of the thermal 
analyzer. These items are automatically entered by the model builder. Item 1 describes the 
accuracy between iterations. The thermal analyzer uses an iterative method to arrive at a 
solution. As soon as the variation of every temperature computed does not differ from its 
previously computed value by the error criterion specified here, computation stops. There 
is a trade off between accuracy and computer time involved in selecting this number. Item 
2 is a damping value which is used to prevent temperature oscillations between iterations. 
Item 3 is the maximum number of iterations to be performed by the computer. This 
number is specified to prevent long computer runs in the event of faulty data. Item 4 is a 
convergence factor. Excessive run times may be induced by the damping factors if there is 


a slow convergence of either overdamped or underdamped temperature values. If it has 
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been determined that this is occurring, the convergence factor will increase the convergence 
due to damping in order to reduce computer time. Item 5 specifies the initial temperature 
estimate for use in the iteration process. 

Input data set 1 contains the temperature dependent coefficients for as many as 15 
curves. These curves and their structure have been previously described in the microcircuit 
modelling section. 

Input data set 2 contains up to 50 constant temperatures. 

Input data set 3 involves heaters for fast warm up and is not used in this model. 

Input data set 4 contains the information for the construction of the n node equations. 
For each node there will be two lines of data. The odd numbered lines specify which 
nodes are connected to the node in question and the mode of heat transfer. Each mode of 
heat transfer is assigned a numerical value from 1 to 10. Our model builder is concerned 
only with conduction which has been designated by the number one (1). For example, if 
node 25 was connected to node 2 by conduction, in line 3 of data set 4 an entry of 251 
would appear. The even numbered lines specify the internodal conductances. For the 
model builder this line will contain the curve number specified in data set 1 corresponding 
to the connecting nodes. 

Input data set 5 specifies the use of unique exponents and is not used in this model. 

Input data set 6 specifies the initial temperature guesses corresponding to the number 
of nodes receiving secondary heat and is not used in this model. All initial temperature 
guesses are equal and are set in line 5. 

Input data set 7 specifies the number of temperature dependent heat input curves and is 
not used in this model. 
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C. FEATURES OF MICCKT 


The MICCKT program will generate a file which can be read by the thermal analysis 
software in order to provide the thermal analysis. The program is user friendly and will 
save numerous man hours in the analysis process. Through the use of MICCKT reaching 
the 750 node capability of the analyzer is now possible with a minimum of effort. 

The program is interactive and will provide prompts for the required inputs. After the 
prompt for a title for the output file there will be a list of options to produce the number of 
nodes for analysis. These options are: 

1) 16 by 16, (640 nodes) 

2) 12 by 24, (720 nodes) 

3) 10 by 30, (750 nodes) 

4) 8 by 32, (640 nodes) 

5) 12 by 18, (540 nodes) 
The options were selected to provide the user with several aspect ratios to approximate the 
microcircuit structure. Option 1 provides the user with a 1:1 aspect ratio while options 2 
through 5 provide the aspect ratios of 1:2, 1:3, 1:4, and 2:3 respectively. 

After selecting the number of nodes for analysis the user must specify the use of either 
English or SI units. Once this selection is made all subsequent inputs will conform to the 
chosen units. The next several inputs are concerned with the physical size of the chip and 
with initial conditions. The user will input the length, width, and overall thickness of the 
chip, the initial temperatures of the lid and sink, and the air gap distance between the chip 
and the lid. As a reminder, the user is limited to an air gap thickness of 0.254 centimeters 
(0.10 inches) maximum. 

It is at this point that any heat inputs to the chip may be entered. The user may select 


from four options. The first option is total rate of heat applied to the chip. The total rate of 
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heat will be divided by the number of nodes in the top layer and evenly distributed to each 
node as an average input. The second option is average rate of heat per unit area, this will 
be either in Btu/hr-in2 or Watts/cm2. This average rate of heat per area will be multiplied by 
the area of the chip and evenly distributed to each node. The third option allows the user to 
select specific nodes for individual heat rate inputs. The fourth option is for no external 
heat input to the chip. 

Once these selections are made a summary of the user’s inputs is displayed for 
verification. The user can review these inputs and make necessary changes before 
continuing. If all entries are satisfactory the run will continue and produce an output file. 

This output file will be in the the proper format for use in the thermal analyzer. All 
data input lines and data input sets will be automatically generated. In particular, for input 
data set 4, each node will be connected to its neighboring nodes by the mode of heat 


transfer and appropiate temperature dependent curve. 
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APPENDIX A 


THIS PROGRAM 15 DESIGNED TO PROVIDE AN OUTPUT FILE WHICH WILL BE 
READ BY THE THERHAL ANALYSIS SOFTUARE. IT 1S A MODEL BUILDER FOR A 
MICROCIRCUIT STRUCTURE OF PURE SILICON AND WILL TAKE INTO 
CONSIDERATION THE VARIATIONS It THERMAL CONDUCTIVITY EXHIBITED BY 
SILICON. THE STRUCTURE IS DIVIDED INTO FOUR LAYERS WITH THE TOP THO 
LAYERS CONTAINING FOUR TINES AS MANY NODES AS THE BOTTOM THO LAYERS. 
THE STRUCTURE IS COVERED HIITH A LID AND THE GAP BETHEEN THE LID AND 
CHIP IS AIR FILLED. 

COMMENTS ARE ENTERED THROUGHOUT THIS PROGRAM Itt ORDER TO EXPLAIN 
THE DEFINITION OF SYMBOLS USED OR TO EXPLAIN THE LOGIC INVOLVED. 
THIS PROGRAM HIAS DEVELOPED AS PART OF A THESIS BY LCDR JOSEPH A. 
WILLHELM, “COMPUTER AIDEN THERMAL ANALYSIS OF A WICROCIRCUILT 


ei here sees? DEC. 1990. FURTHER INFORMATION CAN BE OBTAINED IN THIS 
HESIS. 


CHARACTER ® 79 TITLE 
CHARACTER ® 1 RESP, YES, ANS, CHANS 


THE TITLE ALLOWS FOR UP TO 79 ALPHANUNERIC CHARACTERS TO BE 
ENTERED TO DESCRIBE THE OUTPUT FILE. RESP IS THE RESPONSE TO HAKING 
ANY CHANGES TO THE SUNNARY bIHEN IT APPEARS. YES IS THE RESPONSE TO 
DIFFERENT QUESTIONS AS THEY OCCUR, ANY OTHER RESPONSE WILL RESULT IN 
CONTINUATION. ANS IS THE RESPONSE TO MAKING ANY CHANGES TO THE HEAT 
INPUTS. CHAHS IS THE RESPONSE TO ALL CHANGES MADE IN THE HEAT INPUTS. 


INTEGER R, COUNT, CNTR, NODE, I. J, XROW, YROW, SELECT, NDE, 
+CHNODE, LVL, COUNTER, WX, ty, HR, NL, CTR, KNTR, 
+ALPHA, BETA, GAMA, DELTA, NP, Wl 


INTEGER CONST, HIR, FYP, SECHT, COEFCV, FASTHE., WICV, 
4CONMIROLT, CONTROL2, COUTROL3, NODES, KINAX, WAX, XP, SHT, EFCV, 
+TCV, STHT, MAXITS, Dit, UK, UNIT, LID, SKN 


THE FOLLOWING ARE A LIST OF THE ABBREVIATIONS USED Itt THE INTEGER 


LISTING: 
COUNT, CHTR, COUNTER, CTR, AND KNTR = COUNTERS 
NODE, NDE = NODE NUMBERS OF THE STRUCTURE 
I, J, AND NH = USED IN DO LOOPS, COUNTERS 
XROMN = NUMBER OF NODES IN THE X DIRECTION 
YROW = NUMBER OF MODES IN THE Y DIRECTION 
NX = XROW7Z2, NUMBER OF NODES It THE X DIRECTION OF THE 
BOTTOM THO LAYERS 
NY = YROM72, NUMBER OF NODES IN THE Y DIRECTION OF THE 
BOTTOM THO LAYERS 
NR = NEXT ROW, TO DETERMINE THE NODE ANJACENT TO THE 
NODE UNDER CONSIDERATION, + OR - Y DIRECTION 
NL = NEXT LAYER, TO DETERMINE THE NODE ADJACENT TO THE 
NODE UNDER CONSIDERATION, ABOVE OR BELO! 
R = TOTAL WUNBER OF HODES SELECTED 
HP = NEXT POINT, COUNTER 
X, Y, AND Z = THREE DINENSTON COORDINATE SYSTEM 
SELECT = NUMBER OF THE OPT10N SELECTED FOR NUMBER OF 
NODES TO BE USED 
CHNODE = CHANGE NODE, TO CHANGE THE NODE SELECTED FOR 
HEAT INPUT 
LVL = LEVEL, IN THIS CASE SET TO G, ALTHOUGH NOT 
SPECIFICALLY USED IN THIS PROGRAM IF FUTURE 
APPLICATIONS REQUIRE A VARIABLE LAYER SELECTION THE 
BASIS FOR THIS PURPOSE ARE INBEDDED It! THIS PROGRAM 
ALPHA, BETA, GANMA, AND DELTA = ALGORITHMS USED TO 
CONNECT THE NODES IN LAYER THO TO THE NODES IN LAYER 
THREE DUE TO THE Gil RATIO BETHEL THE LAYERS 
UNIT = RESPONSE TO SI OR ENGLISH UNITS 
SNK = NUMBER OF SINK NODES 
LIDN = NUMBER OF NODES TO THE LID AND AIR GAP 
ALL OTHER INTEGER ENTRIES INVOLVE NECESSARY INPUTS FOR THE 
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lo) 
HANnANANAAANANAAAAAAAAAANAAAAAAANAAAAAAAAAAAAAAAAAAAAAAAAG 


THERMAL ANALYZER ANID ARE USED HERE FOR ACCOUNTING PURPOSES. 
INTEGER BC(750,800), THPC10), KOUNTERC750), SUBI(9) 


B ARRAY = NODE NUMBER ARRAY 


THP AND SUB1L ARRAYS = TEMPORARY ARRAYS USED FOR COUNTING 
KOUNTER ARRAY = USED AS A COUNTER 


REAL WATTS. TOTHT, HTPAR, AVG, AVG], LENGTH, WIDTH, EPSL, EPSD, 
+EMS, TX, TAY, HWOAP, HRAD, FE, LID, SINK, CHOICE, LINE, CHHATT, 
+THICK, THI, TH2, DELX, DELY, DELZL, DELZ2, RAT, FACTOR, 
+DELZ3, EBELZ4, GAP, TEMP], TENP2, TEMP3, TEMPS, TEMPS, TENP6, 


WATTS = HEAT INPUT TO THE NODES, CAN REPRESENT HEAT INPUT 
IN EITHER WATTS OR BTUCHR 

TOTHT = TOTAL HEAT APPLIED TO THE CHIP 

HTPAR = HEAT PER UNIT AREA APPLIED TO THE CHIP 

AVG = TOTAL HEAT DIVIDED BY THE NUMBER OF NODES IN THE 
FIRST LAYER, PROVIDES EQUAL DISTRIBUTION OF HEAT 

AVG = HEAT PER AREA TINES THE AREA DIVIDED BY THE NUMBER 
OF NODES IN THE FIRST LAYER 

LENGTH = LENGTH OF THE CHIP, X DIRECTION 

WIDTH = WIDTH OF THE CHIP, Y DIRECTION 

EPSL = ENUISSIVITY OF THE LID 

EPSD = EMISSIVITY OF THE CHIP (DIE) 

EMS = EWISSIVITY EQUIVALENCE, 17EPSD + 1ZEPSL - 1 

FE = INVERSE OF ENS, EMISSIVITY FACTOR 

HRAD = RADTATION COUSTANT OF PROPORTIEONALITY FOR ATR 

HGAP = CONDUCTION CONSTANT OF PROPORTIONALILTY FOR AIR 

RAIL = RATIO OF HRAD TO NGAP, CORRECTION FACTOR 

FACTOR = 14RAT, NULTIPLIER TO PUT CORRECTION INTO THE 
THERNAL CONDUCTIVITY CURVE FOR AIR 

TX = STARTING TEMPERATURE TO CALCULATE THERMAL 
CONDUCTIVITY OF AIR TO GENERATE CURVE 1109 

TAV = AVERAGE TEMPERATURE OF THE LID AND IX 

Kl THROUGH K8 = THERMAL CONDUCTIVITY CURVES OF SILICON 
REFER TO THESIS FOR HORE DETAILED EXPLANATION 

KAIR = THERMAL CONDUCTIVITY CURVE OF AIR €1109) 

LID = TEMPERATURE OF THE LID 

SINK = TEMPERATURE OF THE SINK 

CHOICE = CHOICE OF THE NETHOD OF HEAT INPUT 

LENE = LINE NUMBER THAT WILL BE CHANGED AFTER SUMMARY 

CHLATT = CHANGE THE VALUE OF THE HEAT INPUT 

THICK = OVERALL THICKNESS OF THE CHIP 

THL = THICKNESS OF THE TOP LAYER, PRESET AT 0.002 It 
¢€0.00508 Ci) 

TH2 = THICKHESS OF THE SECOND, THIRD, AND FOURTH LAYERS 
(THICK-TH1)73 

DELX = LENGTH DIVIDED BY THE NUMBER OF NODES IN THE X 
DIRECTION TO DETERMINE SEPARATION DISTANCE 

DELY = WIDTH DIVIDED BY THE NUNBER OF NODES Ith THE ¥ 
DIRECTION TO DETERMINE SEPARATION DISTARCE 

DELZ1 = NODE SEPARATION IN THE 2 DIRECTION, LEVEL 1 

DELZ2 = NODE SEPARATION 18 THE Z DIRECTION, LEVELS 2-4 

DELZ3 = HODE SEPARATION IN THE Z DIRECTION, BEWICEN 
LAYERS 2 AND 3 

DELZ4 = NODE SEPARATION It THE Z DIRECTION, BETHIEEN 
LAYER @ AND THE SINK 

GAP = DISTANCE BETWEEN THE TOP OF THE CHIP AND THE LID 
PRESET TO A MAXIMUM OF 0.10 IN (0.254 CI) 

TEMP) THROUGH TENP9 = TEMPERATURES SELECTED FROM THE 
SILICON THERMAL CONDUCTIVITY CURVE TO ESTABLISH THE 
CURVES FOR USE It THE THERMAL ANALYZER 


ANSH = RESPONSE TO THE OPTIONS OFFERED FOR NUNBER OF 
NODES TO BE USED 


LIDK = LID TEMPERATURE IN KELVIN 
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Atl QTHFR REAL ENTRIES INVOLVE NECESSARY INPUTS FOR THE 
IHERUAL ANALYZER AND ARE USED HERE FOR ACCOUNTING PURPOSES. 


aaqgagng 


REAL AC750,800), TEMPC1O), LIDTC30), SINKT(4), 
+THCONDC9), CURVEC8,9), CC8), SUBC9), THAIR(9), KARC9) 


A ARRAY = HEAT INPUT ARRAY 

TEMP ANID SUB ARRAYS = TEMPORARY ARRAYS FOR COUNTING 

LIDT ARRAY = INPUTS CONSTANT TEMPERATURE OF THE LID 
FOR THE THERMAL ANALYZER BASED OW REQUIREMENTS 

SINKT ARRAY = IHPUTS CONSTANT TEMPERATURE OF THE SINK 
FOR THE THERMAL ANALYZER BASED ON REQUIREMENTS 

THCOND ARRAY = THERMAL CONDUCIVITY OF SILICON AT THE 
TEMPERATURES SELECTED IN TEMP1 - TEMPS 

C ARRAY = AREA DIVIDED BY THE PATH LENGTH BASED ON INPUTS 

CURVE ARRAY = C ARRAY TINES THE THCOND ARRARY 
THE TEMP1L - TEMPO POINTS IN CONJUNCTION WITH THE 
CURVE ARRAY PROVIDE THE THERMAL ANALYZER WITH THE 
DATA TO GENERATE THE 1101 - 1108 CURVES 

“THAIR ARRAY = EQUATION TO GENERATE THERMAL CONDUCTIVITY 
OF AIR GIVEN A TENPERATURE 

KAR ARRAY = THE K VALUES FOR THE CURVE 1109 


PARAHETER CYES='Y') 


2N ANAAANAANAAAVDAAHHAANHAN 


PRINT *, "ENTER THE WANE OF YOUR FILE.’ 
READ ¥, TITLE 


c 
PRINT *,*YOU MAY SELECT ONE OF THE FOLLOWING OPTIONS.* 
30 PRIUE ¥,'L.) 16 BY 16, YEELDS 640 NODES* 
PRIME *,'2.) L2 BY 26, YIELDS 720 NODES! 
PRINT *,'3.) 10 BY 30, YIELDS 750 NODES* 
PRINT %,'G.) 8 BY 32, YIELDS 640 NODES' 
PRINT x,"5.) 12 BY 18, YIELDS 540 HODES" 
READ x, ANISH 
C 
Cc THESE OPTIONS WILL ALLOW THE USER TO SELECT THE NUMBER OF NODES 
C FOR ANALYSIS. THESE OPTIONS HERE DESIGNED TO ALLO! THE USER TO 
C SELECT FROM A LENGTH TO NIDTH RATIO OF 1:1, 1:2, 1:3, 1:4, OR 2:3. 
Cc 
IF CANSH.NE.1 AND. ANSH.NE.2. AND. ANSH.NE.3.AND.ANSH.HE.G.AND. 
+ANSH.NE.5) THEN 
PRINT *,'ONLY OPTIONS 1.,2.,3..4.,0R 5. ARE ALLOWED.* 
GOTO 30 
ELSEIF CANSH.EQ.1) THEN 
XROH=16 
YROW=16 
CONST=18 
ELSEIF (ANSH.EQ.2) THEN 
XROW=12 
YROW=24 
CONST=20 
ELSEIF (ANSH.EQ.3) THEN 
XROW=10 
YROM= 30 
CousT=20 
ELSEIF (ANSH.EQ.4) THEN 
XROW=8 
YROH= 32 
CONST=18 
ELSEIE CANSH.£Q.5) THEN 
XROME LY! 
-YROU=18 
CONST=15 
ENDIF 
Cc 
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anNAaAA-ANAANg 


120 


aa 


PRINT *, ENTER 1.) FOR SI OR 2.) FOR ENGLISH 
4UNTIS.! 


READ x, UNIT 


THE Mele USED FOR INPUTS IN THIS PROGRAM WILL BE: 
S1; CENTINETUERS, WATTS, DEG CEUTIGRADE 
ENGLISH; INCHES, BTU/HR, DEG FAHRENHEIT 


OTHER PARAMETERS USED If THIS PROGRAM WILL AUTOMATICALLY CONFORM 
TO THE CHOSEN UNITS. 


IF CUNIT.NE.1 AND.UNIT .NE.2) THEN 
PRINT *,'TRY AGAIN® 
GOTO 10 


LSE 
CONTINUE 
ENDIF 


IF (UNIT.EQ.1) THEN 


gerne x, "INPUT CHIP DIMENSIONS IN CENTIMETERS ® 


PRINT *, INPUT CHIP DINENSIONS IN INCHES’ 
ENDIF 


PRINT *%,*ENTER LENGTH: * 
READ *, LENGTH 


PRINT ¥*, "ENTER WIDTH: ' 
READ *, WIDTH 


PRINT X, ‘ENTER THICKNESS: ' 
READ *, THICK 


PRINT %, "ENTER THE THICKNESS OF THE AIR LAYER BETWEEN THE LID* 
PRINT *,'AND THE TOP OF THE DI ee 
READ *, GAP 


IF CUNTT.EQ.1.AND.GAP.GT.0.254) THEN 


en ee LAYER THICKNESS EXCEEDS MAXIMUM ALLOHED (0.254CM)* 
GULU 120 
ELSEIF CUNIT.EQ.2.AND.GAP.GT.0.10) THEN 


ARE aia LAYER THICKNESS EXCEEDS MAXIMUM ALLOWED (0.10 IND* 


ELSE 
CONTINUE 
IDIF 


PRINT *,*ENTER THE TEMPERATURE OF THE LID.' 
READ *, LID 


PRINT *, ENTER THE TEMPERATURE OF THE SINK.! 
READ *, SINK 


PRINT ener THE EMISSIVITY OF THE LID.* 
READ *, EPS 


PRINT *,*ENTER THE EMISSIVITY OF THE DIE.* 
READ x, EPSD 


NX=XROLZ2 
NY=YROMZ2] 

LVL=4% 

NR=XROH 

HL =XROUKYROM 
R=XROMAYROUK] 25 
COUHTER=1 
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ACT, J)=0 
27 CONLTIVE 
26 CONMELNUE 


Cc 
Cc 
Cc INITIALIZES THE MATRICES WITH ZEROES. THIS 15 DONE AS THE SIZE OF 
C THESE MAIRICES CUULD CAUSE EXCESS COMPUTER TINE IF NOT DONE Im THIS 
C FASHION. 
Cc 
Cc 
110 DO 11 I=1,R 
ACI,770)=0.0 
BCI,770)=0 
11 CONTINUE 
e 
Cc IN THE EVENT THE USER WISHES TO MAKE SIGNIFICANT CHANGES TO THE 
C NODE HEAT INPUTS THIS DO LOOP WILL RESET THE NODES AWD HEAT INPUTS 
C TO ZERO. THIS IS DONE IN THE EVENT THAT THE MODE OF APPLICATION HAS 
. CHANGED OR IF ONE OR MORE NODES ARE TO BE ELININATED. 
Cc 
IF CUHIT.EQ.1) THEN 
PRINT X, "HEAT INPUTS HILL BE IN WATTS. * 
TH1L=0.00508 
INTENP=100 
ELSE 
PRINT ®, "HEAT INPUTS WILL BE IN BTU/HR.* 
THI=0.002 
INTEMP=212 
ENDIF 
¢ 
Cc BASED OH THE UNITS SELECTED THIS STATEMENT SERVES AS A PROMPT AID 
C TO SET THE THICKNESS OF THE TOP LAYER AND INITIAL TENPCRATURE FOR 
- THE THERMAL ANALYZER Ith THE PROPER UNITS. 
Cc 
PRIUT *, "HEAT INPUT TO THE MICROCIRCUIT CAN BE ENTERED Jt THES 
CRU XZ PEOLLLOULNG MAMHER. CHOOSE 1.) TO CNTER AS TOTAL HEAT! 
PRINT «, "APPLIED TO THE CHIP, 2.) TO ENTER AS THE AVERAGE HEAT' 
PRINT ¥%,"PER AREA, 3.) TO ENTER NODE BY NODE, OR 4.) FOR NO HEAT®. 
PRINT ¥,'INPUT.! ‘ 
c READ *, CHOICE 


IF (CHOICE.EQ.1) THEN 
PRINT %, "ENTER TOTAL HEAT APPLIED TO THE CHIP* 
READ %, TOTHT 
AVG=TOTHTZHL 
DO 12 N=1,NL 
ACH,770)=AVG 
BCtH,770)=9991 
12 CONTINUE 
ELSEIF (CHOICE.EQ.2) THEN 
PRINT *, ‘ENTER AVERAGE HEAT PER UNIT AREA*® 
READ ¥®, HITPAR 
AVGL=ECHTPARALENGTHRMIDTHIZHL 
DO 13 Nel, NL 
ACH, 770)=AVG1 
BCH,770)=9991 


13 CONTINUE 

- 

Cc THE NUMBER 9991 IS RECOGNIZED BY THE THERMAL ANALYZER TO BE A 
C CONDUCTIVE HEAT INPUT FROM AN EXTERNAL SOURCE. THE NODE NUMBER 770 
c IS ASSIGNED TO TRACK THIS HEAT INPUT. 

Cc 


ELSEIF (CHOICE.EQ.3) THEN 


PRINT *,'YOU HAY SELECT HEAT INPUTS TO THE INDIVIDUAL NODES. 
111 PRINT *,*YOU ARE LIMITED TO A MAXINUM OF 256 INPUTS FOR THE? 
PRINT x, "640 NUDE HODELS, 288 INPUTS FOR THE 720 NODE MODEL, 
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PRINT *,*300 INPUTS FOR THE 750 NODE MODEL, OR 216 INPUTS' 
PRINT X,*FOR THE 540 MODE MODEL. * 


PRINT X, "SELECT THE NUMBER OF HODES THAT HAVE HEAT INpuTSs.* 
READ *%, SELECT 


AQAMMAAAN 


THE REASON FOR THESE OPTIONS IS THAT THEY REPRESENT THE NUMBER OF 
NODES It! THE TOP LAYER. THESE ARE THE ONLY NODES ALLOWED 10 HAVE AN 
EXTERNAL HEAT INPUT BY THIS PROGRAM. 
IF (SELECT.GT.NL) THEN 
PRINT %, "NUMBER OF HODES SELECTED EXCEEDS MAXIMUM ALLOWED. * 
GOTO 111 
ELSE 
CONTINUE 
ENDIF 
Cc 


DO 14 IT=1,SELECT 

PRINT X, ‘SELECT THE NODE TO BE INPUTTED.' 
READ *®, NODE 

PRINT %, ENTER THE HEAT INPUT AT THAT NODE.’ 
READ *, WATTS 

ACHODE,770)=HATTS 

BCHODE,770)=9991 

KOUNTERCI)=I1 


14 CONTINUE 


Cc 
ELSEIF (CHOICE.EQ.4) THEN 
GOTO 17 
ELSEIF (CHOICE.HWE.1.AND.CHOICE.NE.2.AND.CHOICE .WE.3.AND. 
+CHOICE.NE.4) THE? 
PRINT *,°TRY AGATU! 
GOTO 110 
ENDIF 
ze CONTINUE 


112 PRINT *, *SUMMARY OF INPUT BATA.* 
IF CUNET. Ea. 1) THEN 


PRINT TCA, 1Gh2.4)'," 1.) LENGTH OF THE CHIP IS (CM)s ', 
Pi "CA,G12.4)",* 2.) WIDTH OF THE CHIP IS (CH): ‘, 
ma *(CA,G12.4)",* 3.) THICKNESS OF THE CHIP IS (CM): ‘, 
enn *CA,G12.4)',* 4.) THICKNESS OF THE AIR LAYER (CM)1 ‘, 

PRINT "CA,G12.4)%,* 5.) TEMPERATURE OF THE LID IS (DEG C): ', 
spat *CA,G12.4)",* 6.) TEMPERATURE OF THE SINK IS (DEG C): *, 
Zeer *CA,G12.4)',* 7.) THE EMISSIVITY OF THE LID IS: 1 
Pe "CA,G12.4)',* 8.) THE EMISSIVITY OF THE DIE IS: ', 


PRINT *,°9.) THE HEAT INPUTS TO THE NODES ARE: 


IF (CHOICE.EQ.1) THEN 


PRINT 'CA,G12.4,A)",* AVERAGE HEAT APPLIED TO EACH NODE IS: ° 
+AVG,' WATTS' 


ELSCIF (CHOICE.E€Q.2) THEN 


PRINT *CA,G12.4,A)",* HEAT APPLIED TO EACH NODE IS» ‘ 
+AVG1,' WATTS! 


ELSEIF (CHOICE.EQ.3) THE! 
DO 15 T=1,0NL 


IF aeode re THEN 


GoTo 
ELSE 
NDE=I 
PRINT 'CA,I3)", © NODE *, NDE 
PRINT *C€A,G12.4,7)", * WATTS',ACI,770) 


| 








i5 
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ENDIF 
COM TNUE 
ELSEIF CCHOICE.FEQ.4) THEN 
PRINT 4,°NO HEAT INPUTS SERE DESIGNATED. ° 
ENDIF 
ELSE 
PRINT *CA,G12.49)',°* 
+LENGIH 
PRINT "CA,G12.9)', 


LENGTH OF THE CHIP IS CIN): 
WIDTH OF THE CHIP IS CIN)1 
THICKNESS OF THE CHIP IS CIN): 
THICKNESS OF THE AIR LAYER CI)» 


1 
2 
PRINT "CA,G12.4)",* 3 
4 
PRINT 'CA,G12.4)",' 5.) TEMPERATURE OF THE LID IS (DEG F): 
6.) 
7 
8 


* . . . 
— ww —_— — 


PRINT "CA,G12.4)'°,! 


PRINT *€A,G12.4)',! TEMPERATURE OF THE SINK IS (DEG F): 


+SINK 

PRINT "C€A,G12.9)%," 7.) THE EMISSIVITY OF THE LID IS: 
+EPSL 

PRINT '€A,G12.94)',* 8&8.) THE EMISSIVITY OF THE DIE IS: 
+EPSD 


PRINT %,°9.) THE HEAT INPUTS TO THE NODES ARE:'* 


IF (CHOICE.EQ.1) THEN 


PRINT *€A,G12.4,A)",* AVERAGE HEAT APPLIED TG EACH NODE IS: 
+AVG,' BTU/HR' 


ELSEIF (CHOICE.EQ.2) THE 


| 
PRINT "CA,G12.4,A)"," HEAT APPLIED TO EACH NODE IS: 
+AVG1,' BIU/HR' 


ELSEIF (CHOICE.EQ.3) THE! 
DO 22 I=1,Hl 


IF (BCI,770).€Q.0) THEN 
GOTO 22 


»', © NODE ', NDE 
2.4.7)", * BTUZHR',ACI,770) 


CONTINUE 
ELSEIF (CHUICE.EQ.4) THEN 


PRINT ®,'NO HEAT INPUTS HERE DESIGNATED. * 
ENDIF 


ENDIF 


PRINT %,*DO YOU MISH TO MAKE ANY CHANGES, Y OR NT 
READ (*%,5) RESP 


FORMAT CAL) 


IF CRESP.EQ.YES) THEN 


PRINT ¥, "WHICH LINE ITEM DO YOU BIISH TO CHANGE?’ 
READ *®, LINE 


LS 
GOTO 113 
ENDIF 


IF CLINE.EQ.1) THEN 


PRINT *, ‘ENTER THE LENGTH OF THE CHIP.' 
READ x, LENGTH 
GOTO 112 


ELSEIF (LINE.EQ.2) THE 


PRIHT *, ENTER THE WIDTH OF THE CHIP.® 
READ *,MIIDTH 


GOTO 112 


ELSEIF CLINE.EQ.3) THEN 


PRINT «, "ENTER THE THICKNESS OF THE CHIP.' 
READ ®, THICK 


GOTO 112 
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ELSEIF CUIMNE.EQ.4) THE 


PUREE », SPNTER THE THICKNESS OF THE AIR LAYER.* 
READ *, GAP 


GOTO 112 


ELSEIF (LINE.EQ.5) THER 
PRINT x, "ENTER THE LID TEMPERATURE. ' 
READ *, LID 
GOTO 112 


ELSEIF CLINE.EQ.6) THEN 
PRINT *, ‘ENTER THE TEMPERATURE OF THE SINK. 
READ *, SINK 
Goro 112 


ELSEIF CLINE.EQ.7) THEN 
PRINT %, "ENTER THE EMISSIVITY OF THE LID.* 
READ *®, EPSL 
GOTO 112 


ELSEIF (LINE.EQ@.8) THEN 
PRINT %,"ENTER THE EMISSIVITY OF THE DIE.' 
READ *, EPSD 
GoTO 112 


ELSEIF CLINE.EQ.9) THEN 
PRINT *,'DO YOU HISH TO CHANGE YOUR CHOICE OF HEAT INPUT?" 
READ ©%,4) ANS 
FORHAT (Al) 
IF CANS.EQ.YES) THEN 
GOTO 110 


PRINT *, "ENTER NODE “AND HEAT INPUT TO BE CHANGED: * 
READ *, CHMODE, CHHATT 
BCCHNODE, 770929991 
ACCHNODE,770)=CHHATT 
PRINT %, "ALL CHANGES COMPLETED?'* 
READ (%,6) CHANS 
FORHAT CAL) 
IF (CHANS.EQ.YES) THEN 
CHUICE=3.0 
GUTO 112 


LSE 

GOTO 114 
ENDIF 
CHDIF 


ELSEIF CLINE. NE.J.AND.LINE.NE.2.AND.LINE.WE.3.AND.LINE.NE.G.AND. 


4+LINE.NE.S.AND.LINE.NE.6 AND. LINE.NE.7 AND.LINE.NE.8. AND. 
4LINE.NE.9) THEN 


PRINT ¥,'TRY AGAIN’ 


GOTO 112 
LSE 
CONTINUE 
ENDIF 
K1=1101. 
K2=1102. 
K3=1103. 
K4=1104. 
KS=1105. 
K6=1106. 
K7=1107. 
K8=1108. 


Kl - K8 ASSIGH THE CURVES 1101 - 1108 TO SPECIFIC NODES FOR USE 
IH THE THERMAL ANALYZER. 2 


THE FOLLOWING SETS OF ALGORITHMS ESTABLISH THE NODE AND THE 


4) 











C NODES WHICH ARE CONNECTED TO IT. THEY ALSO ESTABLISH THE CURVE TO BE 
z USED IN THE ANALYSIS. 
Cc + 
113 DO 31 Y=1,YROH 
é NO 32 X=1,XROW 
4 I=COUNTER 
IF €X-1.EQ.0.AND.Y-1.EQ.0) THEN 
ACI,1+41)= Kl 
ACI, I+NR)= K2 
BCI, L+1L)= lLOxCI4+1)+1 
‘ BCI, T4hR)= LOKCT4Rd41 
ELSEIF (X#1.GT.XROW. AND. Y-1.EQ.0) THEN 
ACI, 1-1)= K 
ACT, T4+tiR)= K2 
BCI,I-1)= 10xCI-1)41 
é BCI, I+4RI= LOXCI+HR)41 
ELSEIF (Y-1.EQ.0) THEN 
ACI,I-l)= Kl 
ACI,I+])= Kl 
ACI, I+HR)= K2 
BCI,I-L)= 1OKCI-1)+41 
BCI, 1+1)= 1OxCI+1)+1 
é BCI, I4HR)= LOXCE +R) +) 
ELSEIF (X-1.£Q.0.AND.Y+1.GT.YROW) THEN 
ACI,I+1)= Kl 
ACI,I-NR)= K2 
BCI, 1+1)= 10*(C1+1)41 
é BCL, L-URI= LOACI-MRI4L 
ELSEIF (X+t1.GT.XROW. AND. Y+1.GT.YROW) THEN 
ACI,1I-1)= Kl 
ACT, I-NRI= K2 
BCI, i-1)= 10x¢I-1)+41 
- BCI,I-NRI= LOXCI-HR)+1 
ELSEIF (X-1.€Q.0) THEN 
ACI,141)= Kl 
ACI, I-UR)= K2 
ACI, I4NR)= K2 
BCI, I+1)= 1LO¥CI+1)+1 
BCI, T-HIRI= LOXCI-tR)41 
c BCI,I+NR)= LOXCI4HRI41 
ELSEIF (X+1.GT.XROV) THEN 
ACI, I-1)= Kl 
ACI,I-NR)= K2 
ACI, I4NR)= K2 
BCI,I-1l)= 10xCI-1)+41 
BCI,I-NR)= LOXCI-tR)41 
" BCI,I+tR)= LOXCI4HRI+41 
ELSEIF (Y+1.GT.YROMWM) THEN 
ACI,1-1)= Kl 
ACI, I+1)= Kl 
ACI, I-tRI= K2 
BCI,I-1)= 10*CI-1)+1 
BCI, I+1)= 1OxC141)41 
‘ BCI, I-NR)= 1LOXCI-NRJ+41 
ELSE 
ACT,I-1)= Kl 
ACI, [+1)= Kl 
ACI, I-tiR)= K2 
ACI, I4NR)= K2 
BCI,I-1l)= 1LOxCI-1)+41 
BCI,I+1)= LOXCI4+1)41 
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BULL L-MKDT LONC EL -KR 
BCT. Leur)s LONCTonRdey 


ENDIF 


ACT, Teuuy: Ky 
BDci SOMk ds Lowchen yey 


COUNTER=COUNTER+ 


CONTINUE 
CONTINUE 


KNTRENL OL eee 
pO 33 Y21,YRON 
pO 34 X#1,XROW 


T=EKNTR 


IF (X-1.E€Q.0.AND.Y-1.€9.0) TH 
ACI, T¢])= KG i a 

T,T¢uR)d: KS5 

aL4)d= LONCES2 pa] 

/LFURI= LONCT EURO] 


+ 
*1)= LOMCTOLI4L 
“URD= 1LOKCT-URIO] 


1 GT -XRON. AND. Y=1.€0.0) THEN 


+ 

32 
NRd= KS 
2) LOeCT-1)0) 
URD= LOMCT4URd42 


POLAR ADV OSOh TE THEN 

3 . . 
y= KS re 

y= 1ONCT-1961 er 

MRIS LORCTI-HRI+1 


LOMCT¢])¢] 
LORCT-Rd+) 
LORCTONRI 4) 


.GT.XROW) THER 
2 K4 

dz KS 

d= KS 

= JONCT-1)41 
d= 1ORCT-RDO] 
d= LOwCL OUR] 
-EQ.0) THER 

® KS 

)4} 

do] 

MRI] 
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C1,1-NL)= K3 
BCI, I-HL)=10%CI-NLD+1 


KNTR=KNTR+1 


CONTINUE 
CONTINUE 


CNTR=2*NL +1 
NP=HtHL+h 
DO 35:-Y=1,NY 
DO 36 X=1,0X 


I=CHTR 


IF (X-1.€Q.0.AtND.Y-1.EQ.0) THEN 
ACT, 1+1)= KG 
ACT, Tttlxd= K5 
HCL, 1LI1)= LO¥KCT HL th 


BCL, LAX) = LOCI 4X) 41 
ELSEIF (X-1.E&Q.0.AND.Y+1.GT.NY) THEN 

ACT, L413= 4% 

ACT, 1-NX)= K5 

BCE, Ttlr= LoOaRCL+ebyel 

BCI, I-HX)= LOXCI-X)41 
ELSEIF (X+1.GT.NX.AND.Y~-1.EQ.0) THEN 

ACI,I-1)= K4 

ACI, IT4+htX)= K5 

BCI,1I-1)= 10xCI-1)+41 

BCI,I4NxX)= LOXCI+NX) 41 
ELSEIF (X+1.GT.NX.AND.Y4+1.GT.NY) THEN 

ACI,I-1)= K4 

ACI, L-tX)= KS 

BCI,I-1)= 10x*(I-1)+1 

BCI, I-tX)= LOXCI-NX)41 
ELSEIF (X-1.E&Q.0) THE! 

ACI, I+1)= K4 

ACI, I-X)= K5 

ACI, I4NX)= K5 

BCI,I+1)= 10xCI4+1)+41 

BCI, I-NX)= 1LORCI-tX)4] 

BCI, Ltuxd= LOMCI4#Nlx)+} 
ELSEIF (CX+1.GT.NX) THEN 

ACI,I-1)= KG 

ACI, I-tX)= K5 

ACT, Ttuxd= KS 

BCI,1-1)= 1LO*CI-1)41 

BCI, I-HX)= LOKCT-UX) 41 

BCI, T4hxd= LOCI +x) 41 
ELSEIF (Y-1.EQ.0) THEN 
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ACI,I-1)= K& 

ACT, I+1)= K& 

ACL, 1Anx)= K5 

BCI, I-1)= 104CI-~-1)+1 

BCI,T41)= 1OKCI4+1)41 

BCE, L4NX)= LORCTANX 41 
ELSEIF (Y+#1l.GT.NY) THEN 

ACI, 1-L)= 4 

ACI, I141)= KG 

ACI, I-NX)= K5 

BCI,I-1)= 10*CI-1)+41 

BCI, L+l)= lLOxXCI+1l)41 

BCI, I-NX)= LOXCI-NX)41 
ELSE 

ACI, I-1)= K4 

ACI, I4+1)= KG 

ACI, I-NX)= K5 

ACI, IT4HX)= K5 

BCI, J-lL)= 1OKCI-1)41 

BCL,I+1)= 1OxCI+41)41 

BCL, I-NX)= LOXCI-NX)41 

BCI, I4hX)= LOxCI+uxdel 
ENDIF 


ALPHA=-HP+X+(CY—-1 3 XCNRFNX) 
BETA=-tIL+X#+CY-1) XCURFUX) 


GAMHA=-HPEXECY-LIECHRINXIFHR 
DELTA=-HL+#X+CY¥-1 ) XCHRENX) +R 


ACI, I+ALPHA)= K6 
ACI, I+BETA)= K6 
ACT, L4GAMNAD= K6 


ACI, TtpEttad= K6 

ACL, Etnt7a)= K7 

ACLOALPHA, 1) = K6 
ACT#BETA,I)= K6 
ACTIGAMHA, T= K6 
ACT+DELTA,1T)= K6 

BCI, 1+ALPHA)= 1LOXCI+ALPHA)+2 
BCT, L+BETAD= LOXCI+BETAD+1 
BCI, I+GAMMAD= LOXCI+GANMAD+1 
BCL, I1+DELTAD= LOXCI+DELTAD+1 
BCL, T4nl7G)= LOXCI FHL 7G) +1 
BDCI4ALPHA,T)= 1014) 
BCI+BETA,1)= LOXMI+1 
BCI+GANMA, I= 1OXI+) 
BCI+DELTA,I)= 10141 


CNHTR=CHTR+1 


CONTINUE 
CONTINUE 


CTR=2¥YNLANL G41 
DO 37 Y=1,NY 
DO 38 X=1,0X 


I=CTR 


IF (X-1.EQ.0.AND.Y-1.EQ.0) THEN 
1)= K4 


CI,1I+ 
CI,I+HX)= K5 

CI,1+1l)= LOXCI+1)+1 
CL, LT txoy= LOKCI4NX)4)] 


I 
I 
I 
I 
¢ 
I 
I 
I 


m 

c 

"“ 
werrm cory 
te ae onl 


WX)= K5 
L)= 10xCI+1)+1 


1.€Q.0.AND.Y+1.GT.NY) THEN 
= KG 
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BCI, I-NX)= LOXCI-NX)41 


ELSEIF (X+#1.GT.NX.AND.Y-1.E&Q.0) THEN 
ACI,1I-1)= KG 
ACT, T4hX)= K5 
BCT, I-1)}= L0*CI-1)+41 
BCL, TAnx)= OMe T4tixy4) 


ELSEIF (X+1.GT.NX. AUD. Y+1.GT.NY) THEN 
ACI,I-1)= KG 

X= KS 

y= 1O*CI-1)+1 

IXd= LOXCI-NX)41 


~1.EQ.0) THE? 
K4 


wap 

nan 
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= 1OxCI+1)+1 
Y= LOKCI-NX)+4+1 
Y= LOMCT+X)+41 


ELS 


5 

5 

*CI-1)+1 
OxCI-NX)+1 
0 

0 


woeorrrpm 
RAR AR ARH 


KCT4NX) 41 
EQ 
K 
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om 
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+ 
I-NX)+1 
I+hX)+1 
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ENDIF 


ACT, I-NL74)= 
BCI, I-NL74)= 


CTR=CTRtEL 


CONTINUE 
CONTINUE 


K7 
LOXCI-NL74)+41 


THZ=CTHICK-THIIZCLVL-1) 
DELX=LENGTIH/XROW 
DELY=HIDTH/ZYROW 
DELZ1=THl 

DELZ2=TH2 
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DELZ3=(THI+TH2)72 
DELZ4=TH272 


KAIR=1109. 
Cc 
Cc 
C THE FOLLOWING ALGORITHM HILL CALCULATE A CORRECTION FACTOR FOR THE 
C THERMAL CONDUCTIVItY OF AIR TO INCLUDE RADIATION CRFECTIS. THIS I5 
C DONE TO SAVE ON COHPUTER RUN TIME IN THE THERMAL ANALYZER. SEE 
e THESIS FOR GREATER DETAILS. 
Cc 
IF CUNIT.EQ.1) THEN 
EPUB berets? 
LIDK=(579)*(LID-32)4+273.15 
c ENDIF 
EMS=C1/EPSL)+C1/EPSD)~2 
FE=17ENS 
TX=250.0 


DO 200 H=1,9 
TAV=0.5*C(LIDK+TX) 
THAIRCHD=7 .76E-SXTAVt2.87E-3 
HGAP=THAIRCN)/GAP 
HRAD=5 , 669E-8XFEXCLIDK¥R24+TXHH2Z RCLIDK+TX) 
RAT=HRAD/HGAP 
FACTOR=1+RAT 
KARCH) =THAIRCN) XDELXXDELYXFACTORXO.0001/GAP 
TX=TX+4+50 

200 CONTINUE 


Cc 
Cc 
Cc THE FOLLONING GROUP OF ALGORITHMS CONNECT THE TOP LAYER NODES TO 
C THE LID NODES CHODES 751-765). 
C THE THERMAL ANALYZER IS LIMITED TO A MAXIMUM OF 27 CONNECTIONS. 
C THEREFORE BY CONNECTING THE TOP LAYER NODES It! GROUPS OF THENTY THIS 
C LIMITATION CAN BE BYPASSED. THE SAME ANALOGY IS APPLIED TO 
: CONNECTING THE BOTTGM LAYER NODES TO THE SINK NODES (NODES 766-769). 
Cc 
DO 62 I=1,20 
ACI, 751) =KAIR 
BCI,751)=7511 
62 CONTINUE 
Cc . 
DO 63 I=21,40 
ACT,752)=KAIR 
BCI,752)=7521 
63 CONTINUE 
Cc : 
DO 64 1=41,60 
ACI,753)=KAIR 
BCI,753)=7531 
64 CONTINUE 
c 
DO 65 1=61,80 
ACI,754)=KAIR 
BCI,754)=7541 
CONTINUE 
BDO 66 I=81,100 
ACT,755)=KAIR 
BCI,755)=7551 
a CONTINUE 
DO 67 3F=101,120 
ACT,756)=KAIR 
BCT,756)=7561 
of CONTINUE 


DO 68 I=121,140 
ACT,757)=KAIR 
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BC1,757)=7571 
CONTINUE 


CONTINUE 
DMM=2*HL4+NL74+1 


DO 83 I=DNM, DMM+20 


ACT,766)=K7 
BCT,766)=7661 
CONTINUE 


DO 84 IT=DMM+21, DMM+40 


ACL,767)=K7 
BCI,767)=7671 
CONTINUE 


IF CANSH.EQ.1 AND. ANSH.EQ.4) THEN 


LIDN=13 
SHKN SG 

DO 72 I=201,220 
ACT,761))=KAER 
BC1,761)=7611 


CONTINUE 
DO 73 1T=221,240 
ACT, 762)=KAIR 
BC1L,762)=7621 
CONTINUE 
DO 74 IT=241,256 
TCheeee aoe e 
1763)=7631 
continue 


DO 85 I=DMM+61,DMM+60 


ACT,768)=K7 
BCI,768)=7681 
CONTINUE 


DO 86 I=DNM+61, DMM+64 


ACT,769)=K7 
BCI,769)=7691 
CONTINUE 


ELSEIF CANSH.EQ.2) THEN 
LIDN=15 


SUKNEG 
DO 75 I=2 


(I 
continue 
DO 92 1=241,260 
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ACT,763)=KAIR 
BC1,763)=7631 
CONTINUE 


DO 93 1=261,280 
ACT, 764) =KAIR 
BCL,764)=7641 

CONTINUE 


DO 94 I=281,288 
ACT,765)=KAIR 
BCI,765)=7651 

CONTINUE 


DO 8? I= pee DMM+60 
AC(T,768 
BCI, 768)2 7681 
CONTINUE 


DO 88 I=DHM+61,DMM+72 
ACT,769)=K? 
BCI,769)=7691 

CONTINUE 


ELSEIF CANS SW.EQ.3) THEN 
CIDN=15 


1 
BCI,761 
CONTINUE 


DO 78 1[=221,240 
ACT,762)=KAIR 
BCT,762)=7621 

CONTINUE 


DO 79 I=241,260 
CT, 763)=KAIR 
(1,763)=7631 

CONTINUE 

DO 80 1=261,280 
ACT,764)=KAIR 
BC1,7649)=7641 

CONTINUE 


DO 81 1=281,300 
ACI,765)=KAIR 
BCT,765)57651 

CONTINUE 


DO 89 I=DMiH+4G1, DMM+60 
“ACI, 768)=K7 
B(I,768)=7681 

CONTINUE 


DO 90 I=DMt+61,DMM+75 
ACT,769)=K? 
BCI,769)=7691 

CONTINUE 


ELSEIF CANSH.EQ.5) THEN 
LIDN=1) 


HC1,761 
CONTINUE 


DO 91 T=DtM+41, DMM+54 
ACT,768)=K7 
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BC1,768)=7681 
CONTINUE 


ENDIF 


- 


20 ae 


tl 
CONTINUE 


DO 19 T=1,SNKN 
SINKTCI)=SINK 
CONTINUE 


~ 
Lad 
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INITIAL SETTINGS FOR THE THERMAL ANALYZER 


QQAMAOK, 


HTR=0 
EXP=0 
SECHT=0 
COEFCV=9 
FASTHT=0 
HTCV=0 
UNK=0 
CONTROL1=0 
CONTROL 2=0 
CONTROL3=0 
NODES=300 
KMAX=50 
HMAX=6 
XP=] 

SHT=2 
EFCV=15 
STUD eS 
TCV=8 
ACC=.05000 
DAMP=.6667 
MAXTTS=12 
CONVFAC= .8000 


DELYXDELZ1)/DELX 
DELXXDELZ1)/DELY 
DELX*DELY)/“DELZ3 
DELYXDELZ2)/DELX 
DELX¥DELZ2)/DELY 
DELXXDELY)/DELZ2 
GX DELXXDELY)/DELZ2 
DELX*DELY )/DELZ4 


QAAQAANANAND 
Fo ae an he he a 
CONAUNDWN 
Net Ne eee ee ee 
WM a 
Po a ee a ee al 


THE TEMPERATURES SELECTED FROM THE SILICON THERMAL CONDUCTIVITY 
CURVE. THE FIRST GROUP IS IN DEG CENTIGRADE, THE SECOND IN 
FAHRENHEIT. EACH TEMPERATURE GROUP IS FOLLOWED BY THE THERMAL 


CONDUCTIVITY OF SILICON AT THAT TEMPERATURE IN EITHER WATTS/CH-C 
OR BIU/HR-IN-F. 


QAaAAMAAAAANAANA 


—_ 
= 
| | 


T.€Q.1) THEN 
271. 


TEHP6=-228. 
TEMV7=-208., 
TEMPB=-173. 
TENP9=26 .85 
THCOHDC1)=0.25 
THCOUDC2)=4.5 
THCUNDC3)=10.0 


~ 

m 

“s 

= 

1 

4 

N 

wi 

w 
NNNNNN NA 
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704 
705 
Cc 


THCOND(4 
THCONDCS 
THCONDCS 
THCOND(7 
THCONDCS 
THCONDCS 
ELSE 
TEMP L = -456 07 
TEMP 225-4948 .87 
TENP3=-441 .67 
TEHMP4=-427.27 
TEMP 5=-405.67 
TEMP6=-378 .67 
TENP7=-342.67 
TENP8&=-279.67 
TENP9=80.33 
THCOND(1)=1.2 


THCONDC3 
THCOHD(4 
THCONDC5 
THCOND(6 
THCOND(7 
THCOND(S 
THCONDCS 
ENDIF 


DO 23 I[=1,8 
nO 24 J=1,9 


CURVECT, JD=CCI) XTHCOND( J) 
CONTINUE 


CONTINUE 


OPEN A FILE FOR THE OUTPUT. THIS FILE WILL BE IN THE CORRECT 


FORMAT TO BE READ BY THE THERMAL ANALYZER. 


OPEN (3,FILE="MICCKT’, ACORN BORMANLE OEP CUES or REL S 
+RECL=108, STATUS= "NEM! 

WRITE (3,701) TITLE 

FORMAT (€1X,A79) 

HRITE €3,702) R, CONST, HTR, EXP, SECHT, COEFCV, FASTHT, HTCV, 
+UNK, UNIT 

FORMAT (1316) 

WRITE (3,703) CONTROL1, CONTROLZ, CONTROLS 

FORMAT (316) 

rs (3,704) HODES, KMAX, HMNAX, XP, SHT, EFCV, STHT, TCV, UNK, 
+Ut 

FORMAT (€1316) 

WRITE (3,705) ACC, DAMP, MAXITS, CONVFAC, INTEMP 

FORMAT (€2G12.4,112,2612.49) 

no 25 21,8 

HRITE €3,706) TEMPL, CURVECI,1), TEMP2, CURVE(I,2), 
+1LMPS,CURVECI,3), TEMPG, CURVECI,4), TENP5, CURVECI,5), TEMP6, 
+CURVECI,6), TENP7?, CURVECI,7), TEMP&, CURVE(I,8), TEMPS, 
+CURVECTI,9) 

FORMAT (9612.4) 


CONTINUE 


DO 201 N=1,9 


WRITE €3,710) THAIRCN), KARCH) 
FORHAT (9612.4) 
CONTINUE 


| 


WRITE (3,707) (LIDTCH),N=1,12), CSINKTCN),N=1,2) 
707 FORMAT (9612.4) 


DO 50 I=1,R 
COUNT=0 
COUNTER=0 

No 51 J=1,800 
IF (ACL, J). EQ.0) THEN 
GOTO 51 
ELSE 
COUNT =COUNT+1 
TEMPCCOUNT)=ACTI, J) 
ENDIF 
51 CONTINUE 
DU 52 J=1,80 
IF (BCI 
oro 


‘eLSe 
COUNTER=COUNTER# 
THPCCOUNTER)=BC(I,J) 

ENDIF 

52 CONTINUE 


Cc 
WRITE €3,708 
708 FORMAT (15,1 
WRITE €3,709 
709 FORMAT (9612. 
50 CONTINUE 


aes J).EQ.0) THEN 


) 
7, 
) 


Cc 
c 
ENDFILE (3) 
CLOSE (3) 
Cc 
Cc 
STOP 
END 


52 





» CTMPCH),N=1,COUNTER), CSUBICN), N=COUNT+1,9) 
CH),N=1,COUNT), CSUBCN),N=COUNT+1,9) 





APPENDIX B 
The following will demonstrate a typical run using MICCKT. On screen prompts will 
be indicated within brackets ( < > ), responses (Resp) will be directly afterwards. 


Comments will be inserted for clarity. 


< Enter the name of your file > 


Resp: This is an alphanumeric string with a maximum of 76 characters. 


< You may select one of the following options: 
1.) 16 by 16, yields 640 nodes 
2.) 12 by 24, yields 720 nodes 
3.) 10 by 30, yieds 750 nodes 
4.) 8 by 32, yields 640 nodes 
5.) 12 by 18, yields 540 nodes > 
Resp: 3 


< Enter 1.) for SI or 2.) for English. 
Resp: 1 


< Input chip dimensions in centimeters. 
Enter length: > 
Resp: 0.75 


If the English system of measurement was chosen the dimension prompt would be in 


inches. 
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< Enter width: > 
Resp: 0.5 
< Enter thickness: > 
Resp: 0.05 
This will be the overall thickness of the chip. The top layer is automatically set at 
0.002 inches or 0.00508 centimeters. Subsequent layer thicknesses are derived as follows: 
(overall thickness - first layer thickness)/3. 
< Enter the thickness of the air layer between the lid and the top of the die.> 
Resp: 0.2 
The maximum allowable thickness is 0.254 cm (0.10 in). 


< Enter the temperature of the lid: > 
Resp: 40 
Temperatures will be in Centigrade or Fahrenheit depending on the users choice of 


dimensions. 


< Enter the temperature of the sink: > 


Resp: 15 


< Enter the emissivity of the lid:> 


Resp: 0.7 


<Enter the emissivity of the chip:> 


Resp: 0.8 
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< Heat inputs will be in watts.> 
This prompt is a reminder of the units the user is working in. Heat inputs will be in 


Btu/hr if in the English system. 


< Heat input to the microcircuit can be entered in the 
following manner. Choose 1.) to enter as total heat 
applied to the chip, 2.) to enter as the average heat 
per area, 3.) to enter node by node, or 4.) for no heat 
input. 


Resp: 2 


< Enter the average heat per unit area:> 


This will be in either Watts/cm2 or Btu/hr-in2. 


Resp: 5 
< Summary of input data. 
1.) Length of the chip is (cm): 0.7500 
2.) Width of the chip is (cm): 0.5000 
3.) Thickness of the chip is (cm): 0.5000E-01 
4.) Thickness of the air layer (cm): 0.2000 


5.) Temperature of the lid is (deg C): 40.00 
6.) Temperature of the sink is (deg C): 15.00 
7.) The emissivity of the lid is: 0.7000 
8.) The emissivity of the chip is: 0.8000 
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9.) The heat inputs to the node are: 


Heat applied to each node is: 0.62S50E-04 Watts > 
Heat applied was 5 Watts/cm2, this was multiplied by the area of the chip and then 
divided by the number of nodes in the top layer to produce an average heat input to each 
node. If option 1 had been chosen the total heat applied to the chip would be divided by the 


number of nodes in the top layer to get an average input to each node. 


< Do you wish to make any changes, y or n?> 


Resp: y 


< Which line item do you wish to change?> 
Resp: 4 


< Enter the thickness of the air layer between the lid and the top of the die.> 


Resp: 0..25 
< Summary of input data. 
1.) Length of the chip is (cm): 0.7500 
2.) Width of the chip is (cm): 0.5000 
3.) Thickness of the chip is (cm): 0.5000E-01 
4.) Thickness of the air layer (cm): 0.2500 


5.) Temperature of the lid is (deg C): 40.00 
6.) Temperature of the sink is (deg C): 15.00 
7.) The emissivity of the lid is: 0.7000 
8.) The emissivity of the chip is: 0.8000 
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9.) The heat inputs to the node are: 
Heat applied to each node is: 0.6250E-04 Watts > 
Do you wish to make any changes, y or n?> 
Resp: n 
With a no response the program will generate an output file for use with the thermal 


analyzer. 
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